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Abstract. Bulk viscosity of neutron star cores containing hyperons is studied taking into account non-equilibrium 
weak process n + n ^p-|-E~. Rapid growth of the bulk viscosity within the neutron star core associated with 
switching on new reactions (modified Urea process, direct Urea process, hyperon reactions) is analyzed. The 
suppression of the bulk viscosity by superfluidity of baryons is considered and found out to be very important. 
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1. Introduction 

The bulk viscosity of matter in the cores of neutron stars 
has recently attracted great attention in connection with 
damping of neutron star pulsations and gravitational radi- 
ation driven instabilities, particularly — in damping of r- 



modes (e.g., Andersson & Kokkotas 2001). It is well known 
that the bulk viscosity is caused by energy dissipation as- 
sociated with weak-interaction non-equilibrium reactions 
in a pulsating dense matter. The reactions and the bulk 
viscosity itself depend sensitively on the composition of 
matter. 

In the outermost part of the outer neutron star core 
composed mainly of neutrons n with admixture of pro- 
tons p, electrons e, and possibly muons ^ bulk viscosity 
is mainly determined by the reactions of non-equilibrium 
modified Urea process. 



N ■ 



N + p + l + Di, p + N + l^n + N + vi, (1) 



where N stands for a nucleon (n or p), I is an electron 
or a muon, and vi is an associated neutrino. The prob- 
lem of damping neutron star pulsations via modified Urea 
process in npe matter was analyzed long ago by Finzi & 
Wolf ( |1968|) (although the authors did not introduce the 
bulk viscosity explicitly) . The bulk viscosity in npe matter 
was calculated by Sawyer (1989a) and in npefi matter by 
Haensel et al. ( 2000| , hereafter Paper II). 

Deeper in the core, at densities p of a few po (po — 
2.8 X 10^** g cm^'^ is the saturated nuclear matter density). 



direct Urea process may be open (Lattimer et al. 1991 ) 



■ p + I + I'l, p + I 



vi. 



(2) 



It produces the bulk viscosity, which is typically 4-6 or- 
ders of magnitudes higher than that due to modified Urea 
process. T his b ulk viscosity was studied by Haensel & 
Schaeffer (|1992D for npe matter and by Haensel et al. 
( 2000t hereafter Paper I) for npep matter. Note that the 
idea of strong enhancement of vibrational dissipation in 
the neutron stars via a weak process similar to direct Urea 
(beta decay and capture by quasinuceons in npe matter 
containing pion condensate) was put forward by Wang & 
Lu ( |1984D . All the studies cited above are focused on not 
too young neutron stars which are fully transparent to 
neutrinos. We will also restrict ourselves to this case. 

At about the same densities, hyperons may appear in 
the neutron star cores (first of all, and A hyperons, 
and then S", S~, S+). To be specific, we will mainly con- 
sider S~ and A hyperons. Once appeared, the hyperons 
may also initiate their own direct Urea processes (Prakash 
et al. 1992) giving additional contribution to the bulk vis- 
cosity, nearly as high as that due to nucleon direct Urea 
process (||). However, direct non-leptonic hyperon colli- 
sions which go via weak interaction (with strangeness non- 
conservation) such as 



n - 
A- 



n ^ 



P 



N, 



(3) 

(4) 
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are much more efficient. They may increase the bulk vis- 
cosity by several orders of magnitude above the direct- 
Urea level. The effe ct wa s analy z ed by L anger & Cameron 
( |1969| ) and Jones (|l97l| , ^00l4 |2001bD . Analogous effect 
in quark matter was studied by several authors. The en- 
hancement of the vibrational dissipation via non-lepton 
strangeness-changing quark collisions was considered by 
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Wang & Lu (1984). The appropriate bulk viscosity was 
calc ulate d by Sawyer ( |1989bD Madsen ( [1992D , Goyal et 
al. (|199|), and Dai & Lu ( |1996D . 

Calculation of the bulk viscosity limited by non- 
leptonic processes in hyperon matter is a complicated 
problem. There are a number of processes of comparable 
efficiency. The matrix elements can easily be calculated 
in the approximation of bare particles and exact SU(3) 
symmetry, and appear to be nonzero for some processes, 
e.g., (^, but are zero for the others, e.g., (^). However, 
experimental data on the lifetime of A in massive hyper- 
nuclei indicate (e.g., Jones 2001b and references therein) 
that process (^) (with = n) is nearly as efHcient as 
"bare-particle" process (||). Calculation of the matrix el- 
ements for "dressed" particles is complicated and model 
dependent; additional complications arise — even in the 
in- vacuum case — due to the SU(3) symmetry breaking 
(Savage & Walden |l997D . 

Another complication is introduced by superfluidity of 
neutron-star matter. It is well known that neutrons, pro- 
tons and other baryons may be superfluid due to attractive 
part of strong baryon-baryon interaction. Superfluidity of 
neutrons and protons has been studied in numerous pa- 
pers (as reviewed, for instance, by Yakovlev et al 



1999 



and Lombardo & Schulze 2001). Hypcrons can also be in 
s uperfl uid state as discussed, e.g., by Balberg & Barnea 
( 1998 ). Critical temperatures Tc of baryon superfluidities 
are very sensitive to the model of strong interaction and to 
many-body theory employed in microscopic calculations. 
Their typical values range from 10* to 10^" K. They are 
density dependent, and they mainly decrease with p at 
densities higher than several p^. 

The effects of superfluidity of nucleons on the bulk vis- 
cosity associated with direct and modified Urea processes 
in npep, matter were considered in Papers I and II. It was 
shown that the superfluidity may drastically reduce the 
bulk viscosity and, hence, the damping of neutron star 
pulsations. 

In this paper we propose a simple solvable model of the 
bulk viscosity in hyperonic matter (Sect. 2) due to process 
(^ and study (Sect. 3) the effects of possible superfluidity 
of n, p, and S~ on this bulk viscosity. In Sect. 4 we discuss 
density and temperature dependence of the bulk viscosity 
in non-supcrfluid and superfluid neutron star cores. 



tage of this model is that it can be solved analytically. We 
will compare it with other models in Sect. 2.4. 

2.2. Matrix element ofn + n^p + S" 

Let us start with the matrix element M in the "bare- 
particle" approximation. The process is described by two 
diagrams with the states of two neutrons interchanged. 
Accordingly, M = M^^^ -I- Af("), and (^ = c = fee = 1) 

A/(") = -A [up-fx{l + Cj5)un'] [u^l\l + C'75)U„] . (5) 

In this case Ui is a standard bispinor, Ui is its Dirac con- 
jugate (i = n, n' , p, S; UiUi = 2mi, where rrii is a bare- 
particle mass), 7^ is a Dirac's gamma- matrix, and 



Gf 

A = -= sin 9c cos 9c 

v2 



(6) 



Furthermore, Gp — 1.436 x 10 ''^ erg cm"^ is the Fermi 
weak coupling constant; 9c is the Cabibbo angle (sin 0c = 
0.231); C = F + D, C = F - D, where D w 0.756 and 
F « 0.477 are the reduced symmetric and antisymmetric 
coupling constants (e.g., Prakash et al. 1992| ). 

Using the standard technique in the limit of non- 
relativistic baryons we sum |Mp over particle spin states 
and obtain 



E 

spins 



|A/|2 = 64|^Px"^'mpmE, x=(l + 3GG')'. (7) 



This expression coincides with that which can be deduced 



from the recent results of Jones (20011:). In the previous 
papers Jones (1971, 2001a) reported analogous expression 
but with x' = (1 — 3GG')^ instead of %. Numerically, re- 
placing minus with plus makes a great difference due to al- 
most total compensation of the terms mx'- X~ 0.001 and 
x' ~ 4.13. Because of the strong compensation we cannot 
rely on the bare-particle approximation. Let us assume 
that a more evolved calculation based on dressed-particle 
technique will lead to the same Eq. (|^) but with the value 
of X renormalized by medium effects. Accordingly we will 
treat x as a free parameter and, to be specific, we will set 
X = 0.1. 



2. Bulk viscosity of non-superfluid matter 

2.1. Model 

Consider non-superfluid hyperonic stellar matter in the 
core of a neutron star pulsating with a typical frequency 
uj ~ 10'^ — lO'' s~^. In the presence of hyperons the contri- 
bution of direct Urea and modified Urea processes, (|^) and 
(|l|), into the bulk viscosity may be neglected. It is sufficient 
to include the non-leptonic weak-interaction processes (^ 
and (H). For the sake of simplicity, let us take into ac- 
count process (^) alone although we assume that matter 
may contain not only but other hyperons. The advan- 



2.3. Non-equilibrium rate 

Due to very frequent interparticle collisions, dense stel- 
lar matter almost instantaneously (on microscopic time 
scales) achieves a quasiequilibrium state with certain tem- 
perature T and chemical potentials fii of various par- 
ticle species i. Relaxation to the full thermodynamic 
( "chemical" ) equilibrium lasts much longer since it realizes 
through much slower weak interaction processes. 

In the case of process (||) the chemical equilibrium im- 
plies 2 fin = fJ-p-^- fJ-'s- In the chemical equilibrium the rates 
[cm~^ s""'^] of the direct and inverse reactions of the pro- 
cess are balanced, F = F. In a pulsating star, the chemical 
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equilibrium is violated (F ^ T) which can be described by 
the lag of instantaneous chemical potentials, 



(8) 



We adopt the standard assumption (e.g., Sawyer 1989a) 
that deviations from the chemical equilibrium are small, 
\'q\ <ti T. If so we can use the linear approximation 



Ar = r - r 



-Xr], 



(9) 



where A determines the bulk viscosity (S ect. 2.4 ). Our 
definition of A is the same as in Sawyer ( 1989a ). Thus 
defined, A is negative. 

Let us calculate the rate F of the direct reaction, nn —f 
of the process. In the non-relativistic approximation 
we have (?i = c = = 1): 



r = 



dpr, 



Ml 



dpp 



2m„(27r)3 2TO„(27r)3 2mp(27r)3 2mE(27r)3 

spins 

X (5(e„ + e',^ -£p- £s) 5{Pn +p'n - Pp - Ps), (10) 

where pi is the particle momentum and £i is its energy. 
The symmetry factor i before summation sign excludes 
double counting of the same collisions of identical neu- 
trons; /i = {1 + exp[(ei — iii)/T]}^^ is the Fermi-Dirac 
function. 

Evaluation of F is standard (e.g. Shapiro & Teukolsky 



1983) and takes advantage of strong degeneracy of react- 



ing particles in neutron star matter. The multidimensional 
integral is decomposed into the energy and angular inte- 
grals. All momenta Pi are placed on the appropriate Fermi 
spheres wherever possible. Introducing the dimensionless 
quantities 



Xi 



1 



we can rewrite the reaction rate as F = F(°^/, with 

4 



n 



-)-oo 



dxi f{xi) 




(11) 



(12) 



where the blocking factors (1 — / (a;)) are transformed into 
the Fermi-Dirac functions f{x) by replacing integration 
variables x — s- — x, and the typical reaction rate F*^°) is 
defined as (in ordinary physical units) 



F 



(0) _ 



(27r)5 h 



G| sin^ 9c cos^ ( 



*2 



m* m*^ pFS T 



2.15 X 10 



38 



mpTTij: Vlfm 



1/3 



T^X ciR-^ s"^ 



(13) 



In this case m* is an effective baryon mass in dense mat- 
ter, pfs is the Fermi momentum of hyperons, ns is 



their number density, Tg = T/(10^K). Note that in Eq. 
(|l3|) we have used the angular integral calculated under 
the assumptions pps + PFp < 2pFn and pps < PFp which 
are usually fulfilled in hyperonic matter [p-pi being Fermi 
momentum of particle species i). 
The integral /, Eq. (|l2|), is: 

e« 



(47r2 + ^2 



(14) 



e« - 1 6 

The ratef = F(o)/ofthe inverse reaction, E p ^ nn, 



is obtained from F, Eqs. ( |10| ) and (12), by replacing ^ 
Then for |^| < 1 



AF = F(") al 



3 ^ 



Finally, from Eqs. 



and (hSh we obtain 



|A| 



(15) 



(16) 



In non-superfluid matter |Ao| = 2Tr'^T^'^^ / {3kBT) . 
2.4. Bulk viscosity 

The bulk viscosity due to the hyperon process (||) is 
calculated in analogy with that due to the modified or 
direct Urea process (Sawyer 1989a; Haensel & Schaeffer 



1992|). The resuh is 
, CW 1 



LUU}, 



|A|B2 l + a2' - \x\B' 
where ni, is the number density of baryons, and 



B 



drj 



C = Ub 



drj 



1 dP 



(17) 



(18) 



dXy,' dub Hb dX^' 

In this case P is the pressure and Xs = n^/ub is the 
fraction of E^ hyperons. The quantities B and C can be 
calculated numerically for a given equation of state. 

The bulk viscosity depends on the frequency w of neu- 
tron star pulsations. Using the results of Sect. 2.3 the dy- 
namical parameter a can be written as 

' TO„ \ rnp 
\m* ) m* 



6.09 



/lOO MeV\ 



rib 
IfaT 



Ifm" 



1/3 



(19) 



For typical values T - 10^ - 10^ 



where oj^ = cj/(10* s^""") 
K, (jj4 ~ 1, rib ^ 1 fm"'^, <C nb, m* ^ O.lnii, B ^ 100 
MeV, X ~ 0.1 we have a 3> 1. Then we may use the 
high-frequency limit in which is independent of B and 
inversely proportional to uj'^: 



Gl 



T-p "T-S C'^XPFS 



3 (2^)3 
X sin^eic cos^ Oc [kBPf 



2.63 X 10 



30 



rn^ 
m„ 



Ifm" 
C 

100 MeV 



1/3 



g cm ^ s ^. 



(20) 
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If, due to interplay of parameters, a ^ 1 one can use more 
general Eq. (^. For instance, we would have a ^ 1 for 
the same parameters as above but at higher temperatures, 
T Z 10^° K (Sect. 4). We could have a ^ 1 even below 
~ 10^° K if the phenomenological constant x is higher 
than the adopted value % = 0.1. 

In the absence of hyperons the bulk viscosity is de- 
termined by direct or modified Urea processes (Sect. 1). 
These processes are much slower than hyperonic ones. 
They can certainly be described in the high-frequency ap- 
proximation in which partial bulk viscosities due to var- 
ious processes are summed together into the total bulk 
viscosity (e.g.. Papers I and II). Thus we will add contri- 
butions from direct and modified Urea processes whenever 
necessary in our numerical examples in Sect. 4. 

Note that all the studies of bulk viscosity of hyperonic 
matter performed so far are approximate. The subject was 
introduced by Langer & Cameron ( 1969 ) who estimated 
dumping of neutron star vibrations but did not calculate 
the bulk viscosity itself. Jones (1971, 2001a) calculated 
eflEective hyperon relaxation times and estimated the 
bulk viscosity but did not evaluate it exactly for any se- 



lected model of dense matter. Recently Jones (2001b) ana- 
lyzed the bulk viscosity of hyperonic matter taking into ac- 
count a number of hyperonic processes but also restricted 
himself to the order-of- magnitude estimates. 

Our approach is also simplified since we take into ac- 
count the only one hyperonic process (^ and neglect the 
others. Even in this case we are forced to introduce the 
phenomenological parameter x (Sect. 2.2) to describe the 
reaction rate. The advantage of our model is that, once 
this parameter is specified, we can easily calculate the bulk 
viscosity (as illustrated in Sect. 4) and introduce the ef- 
fects of superfluidity (Sect. 3, 4). Technically, it would be 
easy to incorporate the contribution of process (^ as well 
as of other hyperonic processes (Sect. 1). However, for any 
new process we need its own phenomenological parame- 
ter (similar to x) which is currently unknown. Generally, 
in the presence of several hyperonic processes, the bulk 
viscosity cannot be described by a simple analytical ex- 
pression analogous to Eq. Nevertheless, in the high- 
frequency limit the contributions from different processes 
are additive and it will be sufficient to add new contribu- 
tions to that given by Eq. (^^. Thus we prefer to use our 
simplified model rather than extend it introducing large 
uncertainties. 



3. Bulk viscosity of superfluid matter 

3.1. Baryon pairing in dense matter 

Now consider the effects of baryon superfluidity on the 
bulk viscosity associated with process (||). According to 
micro scopic theories (reviewe d, e.g ., by Yakovlev et al. 
1999| and Lombardo & Schulze 2001 ) at supranuclear den- 
sities (at which hyperons appear in dense matter) neutrons 
may undergo triplet-state (^P2) Cooper pairing while pro- 
tons may undergo singlet-state (^Sq) pairing. As discussed 



in Sect. 1 microscopic calculations of the nucleon gaps 
(critical temperatures) are very model dependent. Current 
knowledge of hyperon interaction in dense matter is poor 
and therefore microscopic theory of hyperon pairing is 
even much more uncertain. Since the number density of 
hyperons is typically not too large it is possible to expect 
that such a pairing, if available, is produced by singlet- 
state hyperon interaction. Some authors (e.g., Balberg & 
Barnea 1998) calculated singlet-state gaps for A hyperons. 
We assume also singlet-state pairing of S~ hyperons and 
consider the bulk viscosity of matter in which n, p and 
Yi~ may form three superfluids. Since the critical temper- 
atures Tcm Tcp and Tcs are uncertain we will treat these 
temperatures as arbitrary parameters. 

Microscopically, superfluidity introduces a gap S into 
momentum dependence of the baryon energy, e(p). Near 
the Fermi surface (|p — Pf| ^ Pf) we have 



e = ^ - y + t;|(p - pf)^ at p < pp, 
£ = II + J 5"^ + v%{p - pyY at p>pf, 



(21) 



where vp is the Fermi velocity. The gap S is isotropic (inde- 
pendent of orientation of p with respect to the spin quan- 
tization axis) for singlet-state pairing but anisotropic for 
triplet-state pairing. Strict calculation of the bulk viscos- 
ity with anisotropic gap is complicated. We will adopt an 
approximate treatment of triplet-state pairing (with zero 
projection of total angular momentum of Cooper pairs 
onto the spin quantization axis) proposed by Baiko et 
al. ( |2001| ) for calculating diffusive thermal conductivity 
of neutrons. In this approximation the gap is artificially 
considered as isotropic in microscopic calculations but in 
the final expressions it is related to temperature in the 
same way as the minimum value of the anisotropic gap on 
the Fermi surface. 

It is convenient to introduce the dimensionless quanti- 
ties 



T 
TV. 



y = 



S{T) 
T ' 



= sign(x) v^eM- 



y 



(22) 



For the singlet-state pairing (case A in notations of 
Yakovlev et al. 199£) the dependence of ?/ on t can be 
fitted as 



2/A 



= Vl - r 1.456 



0.157 1.764 



while for the triplet-state pairing (case B) 
yB = (^0.7893 + H^") . 



(23) 



(24) 



3.2. Superfluid reduction factors 

We consider the effects of superfluidity on the bulk viscos- 
ity in the same manner as in Papers I and II and omit tech- 
nical details described in these papers. Following Papers 
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I and II we assume that all constituents of matter partic- 
ipate in stellar pulsations with the same macroscopic ve- 
locity (as in the first-sound waves). Then the damping of 
pulsations is described by one coefficient of bulk viscosity 
^. The effects of superfluidity are included by introducing 
superfluid gaps into the reaction rates, T and f , Eq. (p^, 
through the dispersion relations, Eq. (^). These effects 
influence mainly the only parameter A in Eq. (p7|). Quite 
generally, we can write 

X = XoR, (25) 

where Ao refers to non-superfluid matter, Eq. ([I^), and 
i? is a factor which describes the superfluid effects. The 
latter factor depends on the three parameters, R = 
R[yn, Up, 2/s), which are the dimensionless gaps of neu- 
trons, protons, and E~ hyperons. Obviously, i?= 1 if all 
these baryons are normal (yn = 2/p = 2/e = 0). Calculations 
show that one always has i? < 1 in the presence of at least 
one superfluidity. 

Using Eqs. (^7|) and ( p5| ) we can write the hyperon 
bulk viscosity in superfluid matter in the form 



\Xo\B^R 1 



MB 



(26) 



where ao is the non-superfluid value of a given by Eq. (|l9|). 

In the high-frequency limit (Sect. 2.4), which is often 
realized in neutron star matter, we have Ce oc A, i.e.. 



Cs — Co R, 



(27) 



where Co is the bulk viscosity of non-superfluid matter, 
Eq. (po|). Accordingly, superfluidity suppresses the high- 
frequency bulk viscosity. On the contrary, it enhances the 
static {uj = 0) bulk viscosity Cs « 1/A cx 1/R. Moreover, 
superfluidity increases the dynamical factor a and widens 
thus the range of plasma parameters where the bulk vis- 
cosity operates in the high-frequency regime. 

Under our assumptions superfluidity modifies only the 
integral A/ in the factor A given by Eq. ( pj| ) . To generalize 
A/ to the superfluid case it is sufficient to replace Xi — *■ Zi 
in the all functions under the integral in Eq. ([T^). Then 
R can be written as 



R 



AI 



3_d_ 

7r2 




(28) 



in the limit of ^ ^ 0. Here A/q is the value of AI calcu- 
lated for normal matter, Eq. (p^. 

We have composed a code which calculates R numer- 
ically in the presence of all three superfluids. The results 
will be presented in Sect. 4. Here we mention some limiting 
cases in which evaluation of R is simplified. 

3.3. Superfluidity of protons or hyperons 

The cases in which either protons or hyperons are 
superfluid are similar. Let, for example, neutrons and 
E~ be normal while protons undergo ^Sq Cooper pair- 
ing. Accordingly, R = Rp depends on the only parameter 



y — VAp- For a strong superfluidity (r 
y 3> 1) the asymptote is 



Rv 



y 

2 



TT 

"6" 



(29) 



We have calculated Rp in a wide range of y and proposed 
the flt to the numerical data (with the maximum error 
^ 0.5%) which reproduces also the leading term of the 
asymptote, Eq. (E9): 



Rn 



j5/4 ^ ^1/2 



exp (o.5068 - y/oI06^ 



■r 



(30) 



where a^l + 0.3118?/^ and 6 = 1 + 2.556 y^. If hy- 
perons are superfluid instead of protons, the expressions 
for R are the same but y = 2/ae- 

3.4. Superfluidity of protons and S~ hyperons 

If neutrons are normal but protons and E~ hyperons are 
superfluid R = Rpj:{yp,y^) depends on yp = yAp and 
2/E = 2/AS • We have determined the asymptote of Rp^ at 
large yp and j/s. Let Y be the larger gap, y ^maxjj/s, yp}, 
and yQ = mm{yY,,yp}. At Y— yo^VY ^1 the asymptote 
reads 




(31) 



If ?/o ^ then Eq. (^) reproduces the leading term of the 
asymptot e (|2^ ). To prove this one should consider Eq. (31 ) 
at 1 <C VY <C (Y — yo) <C Y and expand the logarithm 



in Eq. (|3T]) in powers of y/Y'^ - y^/Y < f . 

Equation ( |3l| ) becomes invalid ai yo ^ Y . In this case 
Rp's{yp,y^) ~ Rn{Y), where i?„ is described below. 

3.5. Superfluidity of neutrons 

Now let neutrons be superfluid while protons and E~ hy- 
perons not. For a strong superfluidity (r = T/Tcn ^ 1, 
y = J/Bn > 1) we get 



Rn = ^^e ^, 7 



dq dq' {q — q 
e 9'^ - e 



1.413. (32) 



We have calculated _R„ numerically in a wide range of y 
and proposed the fit (with the maximum error ~ 0.2%): 



Rn 



(o.6192 + v/0.38082 -f 0.1561; 
X exp (o.7756 - ^/0T756^^^^ + 0.18766 y^ 
X exp (^1.7755 - yLTTSPTv) • (33) 
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nN<-^pN/^ 
nN^pNe 



J I I I I I I I I I I I L 



0.6 



Fig. 1. Density dependence of partial bulk viscosities as- 
sociated with various processes (indicated near the curves) 
at T = 10^ K and uj = 10"' s^'^ in non-superfluid matter. 
Dotted and dashed lines refer to Urea processes involv- 
ing electrons and muons, respectively; dot-and-dashed line 
refers to hyperon process (||). Thick solid line is the total 
bulk viscosity. 



4. Results and discussion 

4.1. Non-superfluid matter 

For illustration, we use the equation of state of matter 



in the neutron star core, proposed by Glendenning (1985) 
in the frame of relativistic mean field theory. Specifically, 
we adopt case 3 considered by Glendenning in which the 
appearance of n, p, e, /i, S~, and A is allowed. (Note a 
misprint: numerical values of the parameters b and c of the 
Glendenning (1985) model should be replaced as 6 ^ 6/3 
and c c/4). In this model, muons appear at the baryon 
number density nt = 0.110 fm^'^ (at p — 1.86 x lO^** g 
cm~^); A hyperons appear at ni, = 0.310 fm^"^ {p — 5.51 x 
10^'* g cm~^); and hyperons appear at Tif, — 0.319 
fm^"^ {p — 5.69 X lO'^'* g cm^"^). The density dependence 
of the fractio ns of various particles is shown in Fig. 9 of 
Glendenning (|l985| ). Let us remind that saturation density 
of nuclear matter po w 2.8 x 10'^'' g cm^"^ corresponds to 
Ubo ~ 0.16 fm~^. 

Figure 1 shows the partial bulk viscosities and the total 
bulk viscosity versus n;, at T = 10^ K for stellar vibration 
frequency uj = 10^ s~^. One can see three density intervals 
where the bulk viscosity is drastically different. 

At low densities, ni, < 0.227 fm^"^, the bulk viscos- 
ity is determined by modified Urea processes (Paper II). 
For n(,< 0.110 fm^"^ it is produced by neutron and proton 
branches of Urea process involving electrons (processes (|]) 
with N = n OT p and with ^ e). At higher nt muons are 



created and muonic modified Urea processes (P (again 
with N = n or p but now with / = p) introduce comparable 
contribution. Note that Eq. (29) of Paper II for the an- 
gular integral ApiQ of the proton branch of modified Urea 
process is actually valid at not too high densities, as long 
as ppn > '^PFp—PFi- For higher densities, it is replaced with 



(4^ 
4p2 



Fp 



3 

PFr, 



1 

PFp 



(34) 



which was neglected in Paper II (this replacement has no 
noticeable effect on the values of bulk viscosity). 

At intermediate densities (0.227 fm~'^ < ni, < 0.319 
fm^'^) the main contribution into the bulk viscosity comes 
from direct Urea processes (Paper I). As loiig as rif, < 
0.293 fm^'^ the only one direct Urea process (H) operates 
with I — e while at higher Ub the other one with / = p, 
is switched on; it makes comparable contribution. We see 
that direct Urea processes at intermediate densities am- 
plify the bulk viscosity by more than five orders of mag- 
nitude as compared to the low-density case. 

Finally, at high densities {ni, > 0.319 fm~^), according 
to the results of Sect. 2, the bulk viscosity increases fur- 
ther by about four orders of magnitude under the action 
of non-leptonic process (^ involving hyperons. These 
values of the bulk viscosity are in q ualitative agreement 
with those reported by Jones ( 2001b ). If our model of bulk 
viscosity were more developed and incorporated the con- 
tributions of process (jj) involving A hyperons then the 
high-density regime would start to operate at somewhat 
earlier density, at the A hyperon threshold, rif, = 0.310 
fm^'^. The associated bulk viscosity is expected to be of 
nearly the same order of magnitude as produced by S" 
hyperons (Jones 2001b). Actually, in the presence of hy- 
perons, some contribution into the bulk viscosity comes 
from modified and direct Urea processes involving hyper- 
ons (e.g., Prakash et al. |1992| , Yakovlev et al. |200l| ). This 
contribution is not shown in Fig. 1. It is expected to be 
smaller than the contributions from nucleon modified and 
direct Urea processes (|l|) and (||) displayed in the figure. 

Figure 1 refers to one value of temperature, T = 10^ 
K, and one value of the vibration frequency, uj — 10"^ s^^. 
Nevertheless one can easily rescale ( to other T and cu in 
non-superfiuid matter in the high-frequency regime. For 
the modified Urea processes (M), direct Urea processes 
(D), and hyperonic process (E) we obtain the estimates: 



Cm 



5 X IQi^Tl 



Cd 



5 X 10^4 Ti 



10 



30 



^ gcm-is-1. 



(35) 



The difference in magnitudes and temperature depen- 
dence of C comes evidently from the difference of corre- 
sponding reaction rates. It can be explained by different 
momentum space restrictions (different numbers of parti- 
cles, absence or pres ence o f neutrinos) in these reactions 
(e.g., Yakovlev et al. ^OOlj ). 
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IS Tshear ~ 10r| yrs. It is shown in Fig. 2 by the dotted 
Une. This damping dominates at low T while the damping 
by bulk viscosity dominates at higher T. The total vis- 
cous damping time t (t~^ t;^^^^ + t~^^^^) is displayed 
in Fig. 2 by the solid lines (for the three high-frequency 
bulk- viscosity damping regimes). One can easily show that 
damping by bulk viscosity associated with modified Urea 
processes dominates at T 10^ ujV'^ K. For direct Urea 



-,8, ,1/3 



K, and for hyper- 



processes it dominates at T <; 10 
onic processes at T 3 x 10^ 

Finally, let us mention the validity of high-frequency 
bulk viscosity regime. As follows from Eq. (|l7| ) it is valid 
as long as a J; 1, i.e., lo <; Wc, where the threshold fre- 
quency uJc ~ \\\B/ni,. From Eq. (p^ ) for the hyperon 
bulk viscosity we have oS^ ~ 500 xTg s^^. Using the re- 
sults of Papers I and II we obtain uj^ ~ 5 x lO^^Tg s~^ 
and ~ 5 X 10~^Tg s~^ for modified and direct Urea 
processes. Therefore, we always have the high-frequency 
regime for modified and direct Urea processes at typi- 
cal temperatures T <, 10^° K and pulsation frequencies 
LO ~ 10'' s"'^. The same is true for hyperon bulk viscosity 
excluding possibly the case of very hot plasma, T ^ IQ^^ 
K. Notice that in the low-frequency (static) limit C, cx 1/|A| 
and the temperature dependence of the bulk viscosity is 
inverted with respect to the high-frequency case. 

4.2. Superfluid reduction 

As discussed in detail in Papers I and II superfluidity of 
nucleons can strongly suppress the bulk viscosity produced 
by direct and modified Urea processes. Now let us use the 
results of Sect. 3 and illustrate superfiuid suppression of 
hyperon bulk viscosity. 

Figure 3 shows this suppression at rib — 0.48 fm~'^ and 
u) = lO"' . We present partial bulk viscosities produced 
by hyperonic processes, as well as by direct and modi- 
fied Urea processes. The straight dot-and-dashed lines are 
the partial bulk viscosities in non-superfluid matter. The 
striking difference of these bulk viscosities is discussed in 
Sect. 4.1. Solid and dashed lines show partial bulk viscosi- 
ties in matter with superfluid protons (lgTcp[K] = 9.3) 
and neutrons (lgTc„ = 8.5). At T ^ 10^° K the high- 
frequency approximation for the hyperon bulk viscosity is 
violated. One can see the tendency of inversion of the tem- 
perature dependence of C at T '--^ lO'^" K associated with 
the transition to the low-frequency regime (Sect. 4.1). At 
T < Tcp superfluidity reduces all partial bulk viscosities, 
dissipation of stellar pulsations (Langer & Cameron 1969, In the temperature range Ten < T < Tcp, where protons 

are superfluid alone, all the three partial bulk viscosities 
are suppressed in about the same manner. This is nat- 
ural (e.g., Yakovlev et al. 



Fig. 2. Schematic representation of temperature depen- 
dence of viscous relaxation time scale r in non-superfluid 
neutron star cores with different compositions of matter 
at stellar vibration frequency u — 10-* s"^ Three dashed 
lines show the relaxation due to high-frequency bulk vis- 
cosity associated either with modified Urea processes, or 
with direct Urea processes, or with hyperonic processes. 
The dotted line presents the relaxation due to shear vis- 
cosity. Solid lines refer to the total (bulk+shear) viscous 
relaxation for the three regimes. 

Now we can estimate viscous dissipation time scales t 
of neutron star vibrations. A standard estimate based on 
hydrodynamic momentum-diffusion equation yields r ^ 
pR? /(,, where 7?, ~ 10 km is a radius of the neutron star 
core, and p is a typical density. For the leading processes 
of three types in the high-frequency regime we have 



10 w| 



yrs, TD 



100 cj| 



0.001 w| 



(36) 



Schematic representation of the temperature dependence 
of these time scales is shown in Fig. 2 by dashed lines. 
Sharp difference of the dissipation time scales comes from 
different magnitudes of bulk viscosities in various pro- 
cesses. In particular, the presence of hyperons in the non- 
superfluid neutron star core results in a very rapid viscous 



Jones 1971, 2001a, 2001b) 



Great difference of possible bulk- viscosity scales is in 
striking contrast with the shear viscosity limited by inter- 
particle collisions. The shear viscosity rj should be rather 
insensitive to composition of matter being of the sa me or - 
der of magnitude as in npe matter (Flowers & Itoh 1979|) , 



i.e., rj ^ lO'^^Tg^^ g cni~^ s~^. It is independent of the 
pulsation frequency u. The damping time of stellar pul- 
sations via shear viscosity in a non-superfluid stellar core 



1999) since the reactions re- 
sponsible for the partial bulk viscosities contain the same 
number of superfluid particles (one proton). Indeed, there 
is one proton in hyperonic reaction (^ and direct Urea 
reaction (^, as well as in the neutron branch N = n 
of modified Urea reactions (|l]). At lower temperatures, 
T < Ten, where neutrons become superfiuid in addition 
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Fig. 3. Temperature dependence of bulk viscosity at 
rib = 0.48 fm~'^ and lo = 10* s^^ in the presence of pro- 
ton superfluidity with IgTcp — 9.3 and neutron superflu- 
idity with IgTcn = 8.5. Dot-and-dashed lines (from up 
to down): partial bulk viscosities due to hyperonic, direct 
Urea and modified Urea processes, respectively, in non- 
superfluid matter. Associated solid and dashed lines: the 
same bulk viscosities in superfluid matter. Vertical dotted 
lines show Ig Tc„ and IgTcp. 



to protons, the suppression is naturally stronger and be- 
comes qualitatively different for different partial bulk vis- 
cosities since the leading reactions involve different num- 
bers of neutrons. Evidently, the suppression is stronger for 
larger number of superfluid particles (as well as for higher 
critical temperatures T^,). 

Figure 4 exhibits the same temperature dependence of 
the partial bulk viscosities, as Fig. 3, but in the presence 
of superfluidity of n, p, and E~ (IgTcp = 8.8, Ig Tcs = 9.8, 
IgTcn = 8.3). One can see that superfluidity of S" hyper- 
ons strongly reduces the partial bulk viscosity associated 
with hyperonic process. As a result, at T ^ 10^ K the total 
bulk viscosity is determined by direct Urea processes, fn 
this regime one should generally take into account the con- 
tribution from direct Urea processes with hyperons (Sects. 
1, 4.f). However, under the conditions displayed in Fig. 4 
this contribution can be neglected. 

Therefore, sufficiently strong superfluidity of baryons 
may reduce the high-frequency bulk viscosity by many or- 
ders of magnitude. This reduction will suppress very ef- 
ficient viscous damping of neutron star pulsations in the 
presence of hyperons (Sect. 4.1). Accordingly, tuning crit- 
ical temperatures Tc of different baryon species one can 
obtain drastically different viscous relaxation times. 
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Fig. 4. Same as in Fig. 2 but in the presence of proton, 
neutron and E~ superfluids with IgTcp ~ 8.8, Ig Tc„ — 8.3 
and Ig Tcs = 9.8. 



Note that relaxation in superfluid neutron star cores 
may also be produced by a specific mechanism of mu- 
tual friction (e.g., Alpar et al. 
^000 , and references therein). 



1984, Lindblom & Mendell 



If the neutron star core is 
composed of n, p, e (and possibly p), this mechanism re- 
quires superfluidity of neutrons and protons, as well as 
rapid stellar rotation. The fact that the (conserved) parti- 
cle currents are, in the case of a mixture of superfluids, not 
simply proportional to the superfluid velocities, implies 
non-dissipative drag (called also entrainment) of protons 
by neutrons. Dissipation (mutual friction) is caused by the 
scattering of electrons (and muons) off the magnetic field 
induced by proton drag within the neutron vortices. The 
relaxation (damping) time associated with mutual fric- 
tion, Tinf, depends on the type of stellar pulsations and 
the physical conditions within the superfluid neutron star 
core, in particular — on the poorly known superfluid drag 
coefficient. Its typical value r varies from ~ 1 s to ~ 10'' s. 
One can expect that similar mechanisms may operate in 
the superfluid hyperon core of a rapidly rotating neutron 
star. If so these mechanisms will produce efficient damp- 
ing of stellar pulsations. Note, however, that theoretical 
description of mutual friction is complicated and contains 
many uncertainties. 



5. Conclusions 

We have proposed a simple solvable model (Sect. 2) of 
the bulk viscosity of hyperonic matter in the neutron star 
cores as produced by process (^ involving E~ hyperons. 
We have analyzed (Sect. 3) the hyperonic bulk viscosity in 
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the presence of superfluids of neutrons, protons, and E~ 
hyperons. We have presented illustrative examples (Sect. 
4) of the bulk viscosity in non-superfluid and superfluid 
neutron star cores using the equation of state of matter 
proposed by Glendenning (1985). In particular, we em- 



phasized the existence of three distinct layers of the core 
(outer, intermediate and inner ones), where the bulk vis- 
cosity in non-superfluid matter is very different (in agree- 
ment with the earlier results of Jones 1971, 2001a, 2001bD . 
This leads to very different viscous damping times of neu- 
tron star vibrations for different neutron star models (the 
presence or absence of hyperons; the presence or absence 
of direct Urea process). If we used another equation of 
state of hyperonic matter the threshold densities ni, of 
the appearance of muons and hyperons, and the fractions 
of various particles would be different but the principal 
conclusions would remain the same. As seen from the re- 
sults of this paper and Papers I and II, the high-frequency 
bulk viscosities in all three layers may be strongly reduced 
by superfluidity of baryons. A strong superfluidity may 
smear out large difference of the bulk viscosities in dif- 
ferent layers. In addition, it relaxes the conditions of the 
high-frequency regime. 

Our consideration of the bulk viscosity in hyperonic 
matter is approximate since we include only one hyper- 
onic process (^) (Sects. 1, 2.1, 2.4) characterized by one 
phenomenological constant x- It would be interesting to 
undertake microscopic calculations of x- Analogous prob- 
lem of quenching the axial-vector constant of weak inter- 
action in dense matter has been considered recently by 
Carter & Prakash ( pOOl ). It would also be important to 
determine analogous constants for other hyperonic reac- 
tions (e.g., for (^)) in the dressed-particle approximation. 
This would allow one to perform accurate microscopic cal- 
culations of the bulk viscosity of hyperonic matter. 

In Sect. 4.1 we have presented simple estimates of typ- 
ical bulk viscosities and associated damping time scales 
of neutron star vibrations in different non-superfluid neu- 
tron star models. Let us stress that the actual decrements 
or increments of neutron star pulsations have to be deter- 
mined numerically by solving an appropriate eigenvalue 
problem taking into account various dissipation and am- 
plification mechanisms (e.g., bulk and shear viscosities; 
mutual friction; gravitational radiation) in all neutron star 
layers, proper boundary conditions, etc. (e.g., Andersson 
& Kokkotas |200l| ). In principle, the vibrational motion 
of various superfluids may be partially decoupled. If so 
our analysis of superfluid suppression of the bulk viscosity 
must be modifled (Sect. 3.2). Nevertheless, the presented 
estimates and the theory of superfluid suppression show 
that one can reach drastically different conclusions on the 
dynamical evolution of neutron star vibrations by adopt- 
ing different equations of state in the neutron star cores 
(with hyperons or without), different superfluid models 
and neutron stars models (different central densities, al- 
lowing or forbidding the appearance of hyperons and/or 
operation of direct Urea processes). We expect that the 
results of this paper combined with the results of Papers 



I and II will be useful one to analyze this wealth of theo- 
retical scenarios. 

Note added at the final submission. After submitting 
this paper to publication we became aware of the paper by 
Lindblom & Owen (2001) devoted to the effects of hyperon 
bulk viscosity on neutron-star r-modes. The authors ana- 
lyzed the bulk viscosity taking into account two hyperonic 
reactions, Eqs. (^) and (^, and superfluidity of S~ and A 
hyperons (but considering non-superfluid nucleons) . Their 
treatment of the bulk viscosity in non-superfluid matter 
is more general than in the present paper since they in- 
clude the reaction (|^). They treat the superfluid effects 
using simplifled reduction factors which is less accurate (a 
comparison of analogous exact and simplified reduction 
factors is discussed, for instance, by Yakovlev et al. 1999). 



The principal conclusions on the main properties of the 
bulk viscosity in hyperonic non-superfluid and superfluid 
matter are the same. 
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